Constraining the Physics of Jet Quenching 
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Hard probes in the context of ultrarelativistic heavy ion colhsions represent a key class of observ- 
ables studied to gain informations about the QCD medium created in such collisions. However, in 
practice the so-called jet tomography has turned out to be more difficult than expected initially. 
One of the major obstacles in extracting reliable tomographic information from the data is that 
neither the parton-medium interaction nor the medium geometry are known with great precision, 
and thus a difference in model assumptions in the hard perturbative Quantum Choromdynamics 
(pQCD) modelling can usually be compensated by a corresponding change of assumptions in the 
soft bulk medium sector and vice versa. The only way to overcome this problem is to study the 
full systematics of combinations of parton-medium interaction and bulk medium evolution models. 
This work presents a meta-analysis summarizing results from a number of such systematical studies 
and discusses in detail how certain data sets provide specific constraints for models. Combining 
all available information, only a small group of models exhibiting certain characteristic features 
consistent with a pQCD picture of parton-medium interaction is found to be viable given the data. 
In this picture, the dominant mechanism is medium-induced radiation combined with a surprisingly 
small component of elastic energy transfer into the medium. 

PACS numbers: 25.75.-q,25.75.Gz 



I. INTRODUCTION 



Hard probes, i.e. high transverse momentum {Pt) pro- 
cesses in the context of heavy ion have long been regarded 
as a useful tool to gain tomographic information about 
the medium produced in ultrarelativistic heavy-ion colli- 
sions [ll-@ . The idea underlying tomography is that high 
Pt Quantum Chromodynamics (QCD) processes can be 
reliably calculated using factorized perturbative QCD 
(pQCD) in vacuum, and since the uncertainty relation 
ensures that the scales of medium physics and hard pro- 
cess are well separated, this implies that high pT par- 
ton production takes place in a heavy-ion (A-A) colli- 
sion with a calculable rate, which is approximately equal 
to the vacuum rate corrected for the number of binary 
nucleon-nucleon collisions in an A-A collision. The effect 
of the medium is then that of a final state interaction, i.e. 
chiefly a modification of the outgoing parton fragmenta- 
tion pattern due to parton-medium interactions. As in 
the case of X-ray tomography, having a known process 
embedded into a medium can then in principle be used 
to image this medium and obtain information about its 
density distribution. 

There are, however, notable differences to X-ray to- 
mography which have prevented this idea from being fully 
utilized so far. First of all, in X-ray tomography the inter- 
action between photon and medium is known to a great 
precision, whereas in high Pt tomography the parton- 
medium interaction depends even qualitatively on the 
a priori unknown degrees of freedom of the medium as 
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probed at the resolution scale of the propagating parton. 
Second, while in X-ray tomography a monochromatic 
photon source with a known position is utilized, in high 
Pt tomography both primary parton momentum and po- 
sition arc only known probabilistically. Third, while the 
medium in X-ray tomography is effectively static at the 
timescalc of photon propagation, in heavy-ion collisions 
it undergoes a violent expansion at speeds of a sizable 
fraction of the speed of light c and hence changes during 
the parton propagation. 

The combination of these effects implies that due to 
the necessary averaging over probabilistically known ini- 
tial variables the sensitivity of observables to the desired 
properties of the medium, i.e. its relevant degrees of free- 
dom and its density evolution, is not as large as one would 
naively expect Q , while at the same time there is an am- 
biguity between model assumption made with regard to 
the bulk density evolution and with regard to the parton- 
medium interaction physics. This is exemplified e.g. by 
the work in — while a detailed comparison of models 
with the data is able to constrain one relevant model 
parameter with good precision assuming an otherwise 
perfect theory, there is no consensus between different 
models as to what this implies for the parton-medium 
interaction physics or the medium density evolution. 

The only way to overcome these obstacles is to study 
the full systematics of both medium and parton-medium 
interaction modelling with a sufficiently large body of 
data. While certain parts of the problem have been in- 
vestigated in various papers, a full overview about what is 
known and what constraints the combined high Pt data 
provide has not been available so far. It is the aim of this 
work to present such a meta-analysis of several system- 
atic studies and to argue that by making use of the large 
body of available data, the features of parton-medium in- 
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tcraction and to a lesser degree also the density evolution 
of the bulk medium can be reasonably well constrained. 

This work is organized as follows: First, an overview 
over experimentally available high Pt observables and 
their characteristics is given along with arguments justi- 
fying the selection of observables discussed in this study. 
In section IIIIl an overview over the different classes of 
assumptions which have to be made when specifying a 
model for hard probes is given, as well with a review of 
what is known about the sensitivity of certain observables 
to these assumptions. Section ITVl then summarized what 
constraints can be derived from selected observables us- 
ing a systematic analysis of a large body of models. Sec- 
tion |V] proceeds by providing supplementary information 
on what is known about several other observables and 
other types of models where currently no systematic re- 
sults arc available and how these results relate to the 
previous section. Finally, a summary of the essentail re- 
sults is provided. 



II. EXPERIMENTAL OBSERVABLES 

High Pt tomography utilizes a large number of observ- 
ables raning from single inclusive hadron suppression to 
fully reconstructed jets. Each of these observables has 
its own advantages and disadvantages, and the following 
list aims to illustrate the gross characteristics, strength 
and weaknesses of each and introduces the relevant ter- 
minology. 



A. Single inclusive hadron suppression 

The least involved observable is the single inclusive 
hadron suppression factor Raa which is defined as the 
yield of high Pt hadrons from an A-A collision normal- 
ized to the yield in p-p collisions at the same energy cor- 
rected for the number of binary collisions. 



Raa{pt-,v) 



TAA{h)daPP/dpTdy' 



(1) 



The normalization with the number of binary collisions 
ensures Raa = 1 in the absence of non-trivial initial 
or final state nuclear effects. The default expectation 
is Raa < 1 in medium since parton-medium interaction 
is expected to lead to a flow of high pT parton momen- 
tum into medium degrees of freedom, thus effectively sup- 
pressing the yield in any given momentum bin. In the 
absence of final state effects, nuclear initial state effects 
can however cause Raa > 1 in some kinematical regions 

Experimentally, Raa can readily be obtained for dif- 
ferent colliding ion species, as a function of collision cen- 
trality, as a function of y/s and, if a reaction plane (or 
Vn event plane where Vn is the nth coefficient in a har- 
monic expansion) is identified, also for given centrality 



as a function of the angle (j) of high Pt hadrons with the 
reaction plane {vn event plane). 

If Raa is obtained as a function of 0, then the spread 
^out = Raa{Q) — Raa{t^I'2') between in plane and out of 
plane emission is an important observable sensitive to the 
medium geometry. Knowledge of i?AA(0) and Raa{t^/'2') 
is exactly equivalent to angular averaged Raa and the 
second harmonic coefficient V2 at high Pt, however a dis- 
cussion in terms of Raa emphasizes that unlike at low 
Pt the asymmetry in momentum space is created by an 
extinction process rather than by pressure gradients as at 
low Pt , therefore we will not discuss high Pt phenomena 
in terms of V2 in the following. 



B. Dihadron correlation suppression 

The direct observables in dihadron correlation mea- 
surements are conditional yields, the so-called yield per 
trigger (YPT) on the near and away side. These depend 
crucially on the momentum range required for a trigger 
hadron. The YPT do not directly reflect any suppression 
of the observed rate of triggered events in A-A collisions 
compared to a p-p reference (which is closely related to 
the nuclear suppression factor Raa)- Thus, for instance 
in an initial state suppression picture where a back-to- 
back parton event is either suppressed or never feels a 
medium, Raa can be at an arbitrarily low value while 
the YPT is unchanged in the medium. 

The near and away side YPT is often binned in mo- 
mentum windows, but sometimes also in terms of the 
fraction of the associate hadron momentum divided by 



trigger hadron momentum, zt = P^'"""^ 



IP!^""^ . Since zt 



has a probabilistic connection to the fractional momen- 
tum z = Ptj^"''^ / p^^*' of a hadron produced from a parton 
with momentum p^"'^* , the away side distribution binned 
in terms of zt is often called D[zt) similar to the frag- 
mentation function. 

When comparing the results from A-A and p-p colli- 
sions, usually the yield ratios Iaa = Ymed/Yvac (where 
Ymed and Yvac are the YPT in medium and vacuum 
respectively) are formed where Iaa ~ 1 indicates the 
absence of a medium modification to the correlation 
strength. Unlike in the case of Raa, there is no strong 
a priori expectation Iaa < 1 in the medium. Various 
biases (discussed below) are potentially capable of gener- 
ating Iaa > 1 for Raa < 1 under the right kinematical 
conditions. 



C. 7-hadron correlations 

If one requires the trigger particle to be a photon and 
correlates any associated away side yield, the resulting 
observable is conceptually more powerful than the di- 
hadron correlations discussed above. In 7-hadron corre- 
lations, photons produced in the hard process itself un- 
dergo neither fragmentation nor final state interactions 
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with the medium, and as a result they carry almost undis- 
torted information about the original kinematics at the 
hard vertex, i.e. zt ~ z and D{zt) ~ D{z) for a photon 
trigger. These relations are exactly true in leading order 
calculations but not at higher orders where hard giuon 
radiation can introduce a momentum imbalance. More- 
over, a photon may also be produced in a parton shower 
rather than directly at the vertex, which further dilutes 
the momentum information. 

Since the dominant leading order production channel is 
the QCD compton scattering qg — > 57 and the anihilation 
process qq — > 177 is combinatorically suppressed, to a 
good approximation the away side parton given a photon 
trigger is a quark. This allows to specifically study the 
quark-medium interaction. 

An obvious experimental disadvantage is that due to 
the smallness of the electromagnetic coupling Ofem rel- 
ative to the strong coupling constant and the need 
to identify photons from the hard vertex out of a large 
number of other photon sources the statistics of 7-hadron 
correlation is usually much worse than for dihadron cor- 
relations. 



D. Fully reconstructed jets 

In vacuum, jets are defined by algorithms combining 
hadrons or calorimeter towers in a certain way. The al- 
gorithms are deisgned to undo the pQCD evolution of 
a parton shower, such that jets on the detector level 
(calorimeter towers) are approximately equal to jets on 
the hadron level (particle tracks) and jets on the par- 
ton level (the output of pQCD calculations). This works 
because in the absence of a medium the flow of energy 
and momentum must always remain inside the jet, and 
allows to study the pQCD dynamics on various levels, 
from global event characteristics such as thrust to differ- 
ential momentum distributions of hadrons inside jets. 

However, this assumption is no longer true in the pres- 
ence of a medium, as the pcrturbatively evolving parton 
shower in general exchanges energy and momentum with 
the non-pcrturbative medium. As a result, there is no 
unambiguous notion what a jet in medium should be: If 
one defines the perturbative part of the shower as jet, 
then the energy of jets will be smaller than the energy 
of the original parton, as part of its energy may be car- 
ried by non-perturbative excitations of the medium. If a 
jet is defined as everything causally correlated with the 
shower-initiating parton, then the jet energy will be in 
general higher than the parton energy, as collisions with 
medium partons may correlate them with the jet, while 
part of their energy is thermal. Finally, if the jet is taken 
to be the flow of energy and momentum of the original 
parton, then the jet energy will be the parton energy, but 
a jet can no longer be defined at the hadron level. 

The application of jet fin ding algorithms in a medium 
faces several complications [Ifll which are related to the 
need to account for fluctuations in the background un- 



derlying the jet. Once jet properties are modified as 
compared to the vacuum case, additional complications 
appear. For instance, in is has been shown that im- 
posing Pt cuts (or to a lesser degree angular cuts) may 
significantly obscure the identification of the original par- 
ton energy scale. 

The most striking manifestation of jet modification in 
a medium is the appearance of monojets which has been 
measured by the ATLAS and CMS collaborations 
and quantified in terms of the dijet asymmetry. The di- 
jct asymmetry is a very complicated observable, as it re- 
quires to account for jet finding and jet energy identifica- 
tion in the presence of a fluctuating background medium 
as well as various biases introduced by the presence of a 
trigger jet in the event. 



III. ASSUMPTIONS UNDERLYING JET 
QUENCHING MODELS 

Usually it is tacidly assumed in the literature that 
whenever a calculation of a particular jet quenching 
model is shown together with experimental data this 
would test the parton-medium interaction model. How- 
ever, this is not the case. 

For the purpose of this discussion, let us define a 
parton-medium interaction model as a way to calculate 
the medium-modified fragmentation function (MMFF) 
A->/.(2,£;,g§|Ti(C),r2(C),...r„(C)), i.e. the distribu- 
tion of hadrons h given a parton i with initial energy 
E and initial virtuality Qq where the hadron energy 
Eh = zE and the parton has traversed a medium along 
the path C where TiiC,) are the medium transport coeffi- 
cients relevant for the process. 

In this formulation, the tomographic information 
about the medium spacetime evolution is now in the C, 
dependence of the Ti , whereas the information about the 
nature of parton-medium interactions is in the Ti de- 
pendence of Di-^h^ where the actual evaluation involves 
suitable integrals of Ti[C,) along the path C,. 

This essential core of the parton medium interaction 
model is not experimentally testable, since we can nei- 
ther prepare initial partons with given {E,Q^), nor can 
we prepare a medium with known and controlled Ti(C,). 
Instead, parton type, initial parton kinematics, medium 
spacetime geometry and initial parton position inside this 
geometry are event by event random variables of which 
we at best know the probability distribution. Any com- 
parison of a parton-medium interaction model thus re- 
quires to embed this model into a larger framework which 
supplies these probability distributions. 

Thus, when we talk about for instance 'testing a ra- 
diative energy loss scenario' there are many assumptions 
made both inside the particular framework of modelling 
radiative energy loss and in the surrounding framework. 
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A. Assumptions in the general framework 

Generically, the framework in which a parton-medium 
interaction model is tested can be characterized as fol- 
lows: The initial hard process is computed using pQCD. 
Apart from the transition from nucleon parton distribu- 
tion functions (PDFs) (e.g. [Illii) to nuclear PDFs 
(e.g. [iB-Ell) to take into account changes in the par- 
tonic structure of the initial state, the hard process itself 
is taken to be identical in vacuum and medium. This 
can be justified by the uncertainty relation comparing 
the inverse hard scale with the inverse temperature of 
the medium to estimate the relevant formation times and 
is at the heart of the idea of tomography: If one could 
not compute the initial parton spectrum before passage 
through the medium reliably, one could not make any 
conclusions with regard to the strength of the final state 
interaction with the medium from the measured spec- 
tra in p-p and A-A collisions and would not have tomo- 
graphic information. 

The hard partons from the primary process arc now 
embedded into a medium at an initial vertex position 
(2^0, 2/0) in the plane transverse to the beam axis with a 
random orientation cf) with respect to the reaction plane 
of the A-A collision. The medium description needs to 
specify Ti{Q for ^ [xq + t cos{<f)),yo + T sin{(l))) where r 
is the evolution proper time. However, in most cases the 
medium evolution model is taken to be a relativistic fluid 
dynamics code with parameters tuned to describe bulk 
medium properties (e.g. [H, H^). In this model 
Ti{T,e, s, . . .) connecting the termodynamical variables 
like temperature T, energy density e or entropy density s 
which appear in the hydrodynamical framework with the 
transport coefficients Ti as needed to obtain the MMFF 
must be specified. 

After propagating the partons through the evolving 
medium, the MMFF can be obtained and convoluted 
with the primary parton spectrum to yield the final 
hadron spectrum. This procedure needs to be averaged 
over all possible initial vertices and parton orientations. 
After the final hadron distribution is corrected for ex- 
perimental cuts, comparison with data can be made. 
Most observables are ratios of in-medium over in- vacuum 
quantities, these also require a baseline computation in 
which the MMFF is replaced by the vacuum fragmenta- 
tion function obtained in the same parton-medium inter- 
action model for the limit of vanishing medium. 

Let us now review a number of model assumptions in 
this chain which are not related to specifics of the parton- 
medium interaction and what is known about their va- 
lidity. 

1. Applicability of LO pQCD for the hard process 

In most models, the initial hard process is computed 
in LO pQCD (supplemented with a phenomenological K- 
factor to account effectively for higher order processes) 



with the underlying assumption that the accuracy of such 
a calculation is sufficient and that the error caused by 
neglecting higher order processes is small as compared to 
other uncertainties in the calculation. 

Since observables such as Raa or Iaa are ratios of 
spectra where any /•sT- factor in the calculation drops out, 
the relevant question is if a LO pQCD calculation can 
describe the shape of hard hadron spectra with sufficient 
accuracy to serve as baseline. If this is the case, then one 
can argue that the hard process itself is described with 
sufficient accuracy and all medium modifications can be 
safely absorbed into the MMFF. In [2l| it has indeed 
been shown that a reasonably good description of high 
Pt hadron spectra can be done with LO-PQCD and a 
A'-factor. 

Note that this question is similar to the validity of 
PYTHIA [2^ for the description of high Pt hadron spec- 
tra. PHYTIA is a LO pQCD framework supplemented 
with the PYSHOW algorithm ^ to model a parton 
shower and the Lund model for hadronization. The latter 
two in combination are nothing but a model to compute 
the fragmentation function and its scale dependence. It 
is known that the simulation of a parton shower cap- 
tures many features of NLO pQCD, but fails in certain 
kinematical situations such as truly hard gluon radiation 
(3-jet events). A similar argument can be made here: 
An LO pQCD baseline with a sufficiently well-computed 
MMFF should be valid to good accuracy except for spe- 
cific event topologies where NLO pQCD is explicitly im- 
portant. 

Energy loss for a primary process computed in next- 
to-lcading order pQCD has been studied e.g. in (25l. [26j. 
however with rather drastic assumptions (see below for 
a discussion) about the medium geometry (hard sphere 
overlap, no transverse expansion) and parton medium 
interaction (mean energy loss only, no fiuctuations) that 
the added value of the NLO baseline is difficult to assess. 



2. The spacetime distribution of the medium 

The spacetime distribution of matter produced in 
heavy-ion collisions is ideally something one would like 
to extract in a model-independent way from the exper- 
imental measurements. However, it is now known that 
this is not possible due to the many averaging steps need 
to compute observable quantities — even simpler inver- 
sion problems turn out to be ill-defined [l^] . 

Thus, any comparison of a model with data requires 
to specify a suitable guess for Ti{(^) a priori and con- 
stitutes a test of (among other things) how well a par- 
ticular combination of parton-medium interaction model 
and medium model describes the data. If one could show 
that the outcome of the computation is largely indepen- 
dent of the medium modelling, then one would not have 
any tomographic sensitivity but could test the physics of 
parton-medium interaction well. If, on the other hand, 
the results of all reasonable parton-medium interaction 
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models would be largely generic given a specific medium 
model, then one could do tomography. In practice, nei- 
ther is the case. 

In [2I] it was demonstrated that one could obtain an 
equally good fit to the Pt dependence of Raa in 200 
AGeV central An- An collisions using the same parton- 
medium interaction models with three medium evolu- 
tion models constrained by bulk data (i.e. incorporating 
longitudinal and transverse matter expansion such that 
bulk low Pt hadron spectra were reproduced correctly). 
However, the extracted medium parameters were a factor 
two different between models. Thus, if one would have 
started with a given model for T!i(e) and simply used the 
energy density as given by the medium model, the three 
scenarios would have resulted in vastly different suppres- 
sion of high Pt hadrons. In the same work, the case of 
an (unphysical) static medium was also analyzed — here 
the medium parameter was a factor 5 different from the 
evolving case. This establishes that simplistic medium 
models (such as e.g. the uniform cylinder decaying after 
a finite time as used in [2^ ) can at best yield an order of 
magnitude estimate of the medium transport coefficients. 

The tomographic sensitivity of high Pt hadrons to 
properties of the evolving medium is further illustrated in 
[20] at the example of the spread between in-plane 
and out of plane emission of high Pt hadrons. It was 
found that the spread shows a sensitivity to the choice 
of the medium evolution model which is of the same 
order of magnitude than the sensitivity to the choice 
of the parton-mcdium interaction model. This means 
that the success or failure of a particular parton-medium 
interaction model to describe data in many cases can 
not be discussed independently from the chosen medium 
model. In general, only combinations of models can be 
tested against the data. Thus, every meaningful test of a 
parton-medium interaction model must involve a system- 
atic study with different well-constrained medium evolu- 
tion models. 



3. Relevance of hydrodynamical timescales 



Hydrodynamical codes (e.g. [19|, |20|, |3l|, |32] ) are a tool 
to evolve a thermalized initial state to a final state. The 
thermalized initial state is usually taken to be present at 
a proper time r^, the equilibration timescale of the sys- 
tem. The endpoint of the hydrodynamical evolution is 
chosen to be an isosurface, often an isothermal at Tp, 
the so-called freeze-out temperature, at which point the 



Coopcr-Frye prescription [33[ is used to convert the hy- 



drodynamical fluid into a distribution of hadrons. 

When embedding hard partons into the background 
of an evolving medium, it is often assumed that and 
Tp are the relevant parameters for the jet as well, but 
there is no compelling reason why this should be the case. 
Hard partons do not probe the degree of isotropization in 
the soft medium or pressure gradients, rather they probe 
the density of coloured scattering centers at a certain 



resolution scale. Thus, a hard parton my well scatter 
from a pre-equilibrium state, or interact with a hadron 
which is free-streaming from a hydro-perspective. 

This problem was investigated in [30[ with the result 
that for the most common class of energy loss models 
(radiative energy loss with a decoherence condition, see 
section fill Bp . the uncertainties are parametrically small. 
In practice, they were found to be 0(15%) in the extrac- 
tion of medium parameters and the spread in i?^^(0), 
i.e. smaller than the variation between different hydro- 
dynamical models 0(100%). 

More advanced hydrodynamical codes do no longer use 
the Cooper- Frye prescription to compute the final hadron 
spectrum, but rather run a hadronic cascade following 
the hydrodynamical decoupling (see e.g. [lO, H^])- Such 
models are superior in describing bulk data, but to date 
there is no consistent scheme how to embed hard partons 
into the hadronic rescattering stage. 



4- Transport coefficients and thermodynamics 

When embedding a parton-medium interaction model 
into a hydrodynamical medium, a model for the relevant 
transport coefficients based on the values of hydrodynam- 
ical parameters needs to be specified. The effect of this 
model choice on the extraction of transport coefficients 
has been studied in [3^. 

For radiative energy loss, the relevant transport coef- 
ficient is q (the virtuality transfer from medium to par- 
ton per unit pathlcngth) which is commonly modelled as 



m = K ■ Q(C)(coshp(C) - sinhp(C) cos «(C)) (2) 

where Q(C) is the local density of scattering centers, 
p(C) is the local transverse fiow rapidity of the medium 
(in the commonly assumed Bjorken-type longitudinal 
flow fields, a hard parton is always longitudinally co- 
moving with the medium at any given rapidity slice and 
never sees a flow difference with the medium) and a is 
the angle between parton trajectory and flow direction. 
K is an overall free parameter to adjust the strength of 
the parton-medium interaction given a certain density of 
scattering centers. 

Common choices for Q are e^^'^,T^ or s. For an ideal 
gas, all these prescriptions are identical and count the 
entropy content of the medium, whereas the second part 
of the equation corrects for the Lorentz contraction of 
the volume induced by the rclativistic flow, and hence 
the apparently higher density of scatterers as seen by the 
propagating parton. However, in the presence of a phase 
transition and a hadronic evolution phase, the three pre- 
scriptions lead to sizable differences in the late stage. 
q ^ significantly emphasizes the hadronic phase as 
compared to g ~ e'^/^ [3a |. In the extraction of q, this 
leads to a factor 2 difference. 
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Based on these results, it appears that the ques- 
tion how the parton-medium interaction is treated in a 
hadronic medium is of some importance at the level of 
quantitative parameter extraction. 



5. Smoothness of the medium 

In recent times, it has become apparent that event by 
event (EbyE) fluctuations in the hydrodynamical initial 
state are an important effect. In other words, it matters 
for observable quantities whether one averages first over 
many initial states and computes the final state by hy- 
dro dynamically evolving an average, smooth state, or if 
one evolves many fluctuating irregular initial geometries 
and averages observables over the resulting final states. 
This has driven the development of EbyE hydrodynami- 
cal codes Hi, Hi-Ii^ 

This development naturally raises the question to what 
degree the modelling of partons propagating through a 
medium is infiuenced by choosing a final state rather 
than an initial state averaging procedure. First studies 
in [U \^ using a non-evolving geometry identified two 
main effects, the non-linear response of Raa to 'holes' 
and 'hotspots' in the density distribution and the cor- 
relation of initial hard parton production vertices with 
hotspots with large differences of factors 2 in the 
extracted transport coefficient between initial state av- 
eraged, smooth medium and final-state averaged EbyE 
fluctuating medium. A subsequent calculation using an 
ideal fluctuating hydrodynamics code [11] showed that 
the two effects almost cancel in practice and that the net 
difference in the transport coefficient between smooth 
and averaged medium is rather 0(20%) H^. A large 
reason for this difference to [4l|, |42| is that the medium 
evolution itself removes initial state fluctuations quickly, 
as they cause large pressure gradients. This is even more 
true in a viscous evolution of initial state fluctuations 
(see e.g. [39j). 

There is. however, a residual effect of fluctuations for 
the normalization of Raa for non-central collisions which 
depends on the size scale of the fluctuations [11]. This 
requires further study with more systematics. 



6. Spatial distribution of hard vertices 

The probability distribution for a hard vertex to be 
found in the transverse plane at some position Tq — 
{xo,yo) is usually taken to be the binary collision pro- 
flle for impact parameter b 



, TA{ro+b/2)TA{ro-b/2) 
P^'^o^yo)- j—^^ , (3) 



where the nuclear thickness function is given in terms of 
the Woods-Saxon nuclear density pAif, z) as 

TA{r) = j dzpA{r,z). (4) 

The relevance of this choice becomes apparent only in 
connection with the hydrodynamical initial state. Con- 
sider the case of a cylinder of constant density with the 
Woods-Saxon radius R as medium model. Embedding 
a binary overlap profile in such a distribution creates 
through the tail of the binary profile a halo of partons 
which never encounter the medium at all. Similarly, 
when a more realistic wounded nucleon (WN), binary 
collision (BC) or color-glass-condensate (CGC) profile is 
chosen for the initial state, this choice infiuences the size 
of the halo of non-interacting partons [s^l • 

As mentioned above, in the case of fiuctuating initial 
conditions, it is particularly important that the same set 
of vertices is used in an event to both generate the hy- 
drodynamical initial state and to sample the hard parton 
production points, otherwise the real correlation between 
hotspots in the initial density and initial hard parton po- 
sition is lost [Hl-li^. 

7. Eikonal parton propagation 

In many implementations of parton-medium interac- 
tion, it is assumed that partons probe the medium on 
straight line trajectories. However, interactions with the 
medium in general transfer energy and momentum and 
hence deflect the trajectory of the parton. Generically, 
one would always expect the hard scale to be larger than 
the temperature, px ^ T, and hence the actual amount 
of deflection should be small for sufficiently hard partons. 

In [4^ the eikonal propagation assumption was tested 
explicitly and found to be good on the 10% level as far 
as the origin of observed partons in the final state is con- 
cerned. It should be noted however that non-eikonal de- 
flection of hard partons has no effect on single hadron 
high Pt observables unless the deflection is strong enough 
to probe appreciably different regions in a hydrodynamics 
or parton kinematics. For instance, in a Bjorken picture, 
the loss of partons produced at midrapidity and being 
scattered to forward rapidity would be compensated by 
a corresponding gain from partons produced at forward 
rapidity being scattered into midrapidity. 

Non-eikonal parton propagation is experimentally ac- 
cessible in dihadron correlations. 



8. Summary 

The following table summarizes the uncertainties asso- 
ciated with various assumptions as estimated above (per- 
centages are always given with reference to the lowest 
value observed in a systematic study). We distinguish 
two cases: 1) the extraction of a transport coefficient 
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Ti like g or e (i.e. a measure of the overall strength of 
parton-medium interaction) and 2) a less averaged ob- 
servable such as the spread S™f between in-plane and 
out of plane emission at high Pt, which is the ratio of 
suppression factors in which many systematic uncertain- 
ties cancel. 





q,e 




references 


LO pQCD 


10% 


<5% 


m 


medium model 


100% 


100% 


[28, 30] 


hydro timescales 


15% 


15% 


[30] 


Ti{T,e,s) model 


100% 


<5% 


[35] 


medium smoothness 


20% 


<5% 


[43j 


eikonal propagation 


<5% 


<5% 


m 



Two observations can be made immediately: First, 
there is a hierarchy in the sensitivity: Not all model 
choices affect the results in the same way. To a first 
approximation, we may thus focus on the quantities to 
which the results are most sensitive. Second, by form- 
ing suitable ratios of more differential observablcs, the 
sensitivity can be focused to one assumption class which 
can then potentially be determined from this particular 
observation. 

How to interpret a 100% sensitivity of q or the spread 
to the choice medium model is in the eye of the beholder. 
If one is interested in extracting transport coefficients or 
the physics of energy loss in a model-independent way 
from data, then these numbers imply that this can not 
be done in a meaningful way. On the other hand, a large 
sensitivity to the hydrodynamical bulk medium model 
means that there is tomographic sensitivity in the ob- 
servable and hence a chance to constrain different bulk 
models through high Pt physics. 

Let us stress again that these estimates assume a con- 
strained model for the medium including a realistic time 
evolution and that the sensitivity to model parameters 
is easily an order of magnitude larger if unconstrained, 
static or otherwise unrealistic medium models are used. 
This may explain why unconstrained models allow for 
spectacular effect sizes which decrease dramatically once 
the calculation is repeated in a constrained framework. 



B. Assumptions about parton-medium interactions 

Let us continue with a discussion of assumptions 
needed to create a model for parton-medium interaction. 



The simplest possibility are elastic QCD interactions 
with medium constituents, in which the recoil energy of 
struck thermal partons is lost from the hard parton |45l - 
[soj . Such models have been suggested as an attempt to 
understand the surprisingly large suppression of heavy 
quarks [sij , but the arguments apply equally well to light 
quarks and gluons. Key quantity determining the mag- 
nitude is the effective mass of the recoiling objects — 
under the assumption that the medium consists of static 
scattering centers, there is no collisional energy transfer. 
Modern formulations of the same idea treat the problem 
in MC codes [13] or by solving coupled Fokker-Planck 
equations numerically [52|. 

The distinctive feature about elastic interaction is that 
they are incoherent. If a parton traverses a medium with 
constant density p and looses the mean energy {A.E) iscatt 
in each scattering processes, then the mean energy loss 
after a length L will be 



(AE) = N,,att{AE)i,,att = (7pL{AE), 



(5) 



(with a being the interaction cross section) as long as 
finite energy corrections are not important, i.e. as long as 
(AE) Iscatt can be taken independent of the initial hard 
parton energy E. 

The mean energy loss per unit pathlength is often cast 
into the form of the transport coefficient e and its vari- 
ance into 62 [53| . 

A different possibility to transport energy away from 
a hard parton is medium- induced radiation [H-Q. The 
underlying idea is that as the hard parton traverses a 
medium, interactions change the transverse momentum 
randomly. The transport coefficient q measures the rate 
dQ^ /dx at which a parton acquires virtuality from the 
medium per unit pathlength by this mechanism. A giuon 
from the virtual cloud surrounding the hard parton can 
decohere from the parent and appear as real radiation if 
the interaction with the medium is sufficient to overcome 
the cloud gluon's virtuality 



54| 



Assuming a parent with energy E and a radiated giuon 
with energy Erad and momentum fc^ transverse to the 
parent where E ^ E^ad ^ (soft and collinear emis- 
sion), the formation time associated with the radiation 
process as given by the uncertainty relation is r ^ E/k'^ 
or, assuming the transverse momentum is acquired by 
interactions with the medium, r ^ E/Q^. Since this vir- 
tuality is picked up as = qr, the two expressions can 
be combined to find r ^ Erad/iqT) which can be solved 
for the typical radiated energy Erad- The process is most 
efficient when as much time for decoherence is available 
as possible, for a medium length L and a parton moving 
with the speed of light this implies t = L, thus 



1. Basic interaction structure 

At the core of any parton-medium interaction model 
is an idea by what processes energy is carried away from 
the hard parton, and commonly the models are labelled 
according to this idea. 



Erad - qL^ r^L^ (6) 

since q in a thermal medium is parametrically ^ T^. 
This different power in the pathlength is often used to 
argue that radiative energy loss dominates over collisional 
energy loss. 
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55| , the Arnold- 
58[, the Guylassy- 
591 1 and the higher 



Common model frameworks for this type of 
parton-medium interaction are the Armesto- Salgado 
Wiedemann (ASW) formalism 
Moore- Yaffe (AMY) formalism 
Levai-Vitev (GLV) formalism 
twist (HT) formalism [i,|5^. 

While early radiative energy loss model used the ap- 
proximation that the parent parton has such large energy 
that the change in its kinematics induced by the radia- 
tion of a soft gluon can be neglected, finite energy and 
length effects for medium-induced radiation have now 
been worked out in a number of models including MC 
frameworks where explicit energy-momentum conserva- 
tion at each vertex can be accounted for [60l - l62| . The 
conclusions of all these works agree qualitatively: While 
dependence of energy loss is observed initially, after 
a short timescale (typically 2-3 fm) kinematical and fi- 
nite length correction effectively lead the system back to 
a linear dependence on pathlength. Thus, for realistic 
kinematics, the dependence of radiative energy loss 
does not lead to dramatic differences as compared with 
an incoherent mechanism. 

There is yet a different class of scenarios based on the 
idea that the AdS/CFT correspondence [63l - l65| can be 
used to describe a QCD medium in the strongly coupled 
limit which yields a parametrically different estimate for 
the energy loss of a hard parton (if everything is treated 
in the strong coupling limit, no jets appear in the model). 
In so-called hybrid models, it is assumed that what is ra- 
diated into the medium still comes from the perturbative 
part of the hard parton wave function, but that the in- 
teraction of the virtual gluons with the medium, and the 
way they are freed, is governed by strong-coupling dy- 
namics |66l468l | . 

Here, what acts on the gluon cloud is a longitudinal 
drag force of order F ^ T^, thus the virtuality scale of 
gluons from the hard parton which is accessible after a 
propagation length i is - (i^ • L)^ - T^L'^. If this 
estimate is inserted into the expression for the formation 
time above, parametrically 



rad 



(7) 



is found. Thus, a strong coupling treatment of the 
interaction of the medium with the virtual gluon cloud 
of a hard parton is has both a different temperature and 
length dependence than a weak coupling treatment. 

A very relevant question is how finite energy correc- 
tions would change the pathlength dependence found in 
this class of models. While there does not seem to be a 
detailed framework in which this question has been an- 
swered, comparison with the pQCD induced radiation 
models suggests that the pathlength dependence may be 
weakened after a finite time by at least one power of L. 

Finally, note that realistic models usually combine dif- 
ferent parton-medium interaction scenarios [s^ |6^, [tO] , 
thus the real question is about the relavtive strength 
rather than about which mechanism is true. 



2. Medium-modified showers vs. energy loss 

As stated above, the general outcome of a 
parton-medium interaction model is the MMFF 
D,^hiz,E,Ql\TiiC),T2iO,---TniO) which describes 
the the QCD evolution to a final state in the presence of 
a medium after a highly virtual parton has been created 
in a hard process. In principle, this evolution includes 
the development of a partonic shower as well as the 
hadronization. In practice, exist ing in-medium shower 
codes like JEWEL [7l|, YaJEM [zl [zl or Q-PYTHIA 
[74| use the formation time for hadrons h with energy 
Eh and mass to/j to argue that r ~ Eh/mf^ corresponds, 
for light hadrons or sufficiently energetic hadrons, to a 
scale outside the medium geometry and hence employ 
vacuum hadronization models. It should be noted that 
this approximation is frequently not justified for baryon 
production or subleading shower fragments. 

However, the majority of parton-medium interaction 
model makes the further approximation that it is suffi- 
cient to consider the energy shift of the leading parton, 
i.e. the MMFF is approximated by 



A^h(2, E, QlmiO) « PiAE, E\T,{0) ® D,^h{z, QI) 

(8) 

where Di^h{z, Qq) vacuum fragmentation func- 

tion and P{AE, E) is the probability distribution for the 
leading parton to experience an energy shift AE before 
fragmentation. 

A priori it doesn't seem obvious that this approxima- 
tion works, since the high virtuality scale Qq after the 
hard process ensures very fast partonic branching even 
before a medium can be formed, thus the object inter- 
acting with the medium is almost never a single on-shell 
parton but rather a partially developed shower, and it is 
far from obvious that the vacuum QCD evolution can be 
factored out and carried out after the medium-induced 
energy loss has been considered. 

However, if one convolutes a vacuum fragmentation 
function with a pQCD parton spectrum at RHIC kine- 
matics, one finds a mean z oi 0.7, i.e. for the observed 
single inclusive hadron spectrum, about 70% of the orig- 
inal energy are in the leading hadron. One can combine 
this with a Lund-type hadronization model |24| to esti- 
mate that this implies that on the partonic side more 
than 70% of the shower energy flow through the leading 
parton, thus the single inclusive hard hadron spectrum is 
dominated by events in which there is little vacuum radi- 
ation and the hadron kinematics is driven by the leading 
parton kinematics. In this particular case, it is to good 
accuracy justified to approximate the situation by lead- 
ing hadron energy loss. Note that this argument does 
not apply to different observables such as dihadron cor- 
relations or jets where the energy loss approximation is 
bound to fail for subleading jet fragments. 

From another perspective, this corresponds to differ- 
ent ways of interating the basic parton-medium inter- 
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action process. In in-medium shower evolution codes 
[tiI - It^ I the parton-medium interaction is embedded via 
a spacetime picture of the evolving shower into the iter- 
ation of basic 1 — >■ 2 splitting processes with explicit en- 
ergy/momentum conservation (inclusing the transverse 
momentum) at each vertex. A similar setup is used in 
the more recent version of the HT formalism [5^ where 
the basic parton-mcdium interaction is part of the (lon- 
gitudinal) DGLAP evolution of the fragmentation func- 
tion. As we will see later, there are reasons to think that 
the model choice here mainly affects the distribution of 
subleading hadrons in the fragmentation function. 

The AMY formalism [s^, [Ell starts with the assump- 
tion of a hard on-shell parton, but evolves the whole 
parton distribution created by medium-induced radiation 
with rate equations. Finally, in older formulations of HT 
[^ or the ASW formalism (54 | the assumption of indepen- 
dent gluon emissions is made which allows to convolute 
the gluon spectrum of a single parton-medium interac- 
tion with the Poisson distribution to compute the total 
energy loss probability density. 



E = 20 GeV 

ii->h* 



There are also hybrid models such as MARTINI [75 1 
in which a vacuum shower is evolved till a medium is 
formed, and the energy loss approxiation is then used for 
each of the shower partons present at this point. 

In formulations in which the parton-medium interac- 
tion is formulated as part of an evolving shower from 
some initial hard scale Qq down to a lower scale Qp, 
there is an additional source of pathlength dependence. 
It can be argued that the uncertainty relation permits 
the medium of length L to modify the shower resulting 
from a parton with energy E only down to a scale Qp 
given by [El, [El] 



Qf = VE/L- 



(9) 



This condition amounts to an additional energy and 
length dependence of the MMFF which can't be easily 
cast into an analytical expression and has been used in 
the resummed HT framework (56j and in the MC code 
YaJEM-D [m. 

As evident from Fig. [U the differences between a 
MMFF computed in a leading parton energy loss picture 
and a MMFF from an in-medium shower evolution can 
be quite substantial. Especially at low z, the in- medium 
shower treatment keeps explicitly track of hadron pro- 
duction from he medium-induced soft gluon radiation 
whereas the leading parton energy loss does not. Based 
on the figure, one may wonder why the energy loss ap- 
proximation works at all. The reason is that the MMFF 
for single inclusive hadron spectra needs to be convoluted 
with the pQCD parton px spectrum, and this convolu- 
tion suppresses the high and low z region. It turns out 
that for RHIC kinematics, the MMFF in this example is 
probed dominantly around z 0.6 where any differences 
are small. This suggests that the suppression of single in- 
clusive hadron production in terms of Raa is particularly 
unsuited to probe details of the MMFF and that more 
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FIG. 1: A comparison of medium-modified fragmentation 
functions of a 20 GeV c?-quark as computed in the energy 
loss formulation (ASW) and an in-medium shower evolution 
(YaJEM with radiative energy loss (RAD) or incoherent drag 
(DRAG)) with the condition that all MMFF lead to the the 
same single hadron suppression [t^ . 



differential correlation observables such as Iaa need to 
be studied to access this information. 



3. Mean energy loss vs. fluctuations 

If one assumes that the whole parton spectrum is 
shifted by a mean energy loss {{AE)), one can describe 
the hadron suppression observed at RHIC with a value 
of about {{AE)) ^ 5 GeV [t^. However, such an exer- 
cise is not particularly meaningful, as it extracts a quan- 
tity both averaged over geometry (i.e. the tomographi- 
cal information) and dynamical fluctuations given a path 
in the geometry (i.e. the microscopical dynamics of the 
parton-mcdium interaction), which is in addition biased 
by the fact that a hadron has been observed at high Pt, 
whereas there is no guarantee that all hard partonic scat- 
terings lead to observed hard hadrons. 

In some computations, the geometrical averaging has 
been undone but the dynamical fluctuations are averaged 
over. This corresponds to a fixed value of energy loss 
given a path, but since partons travel on different paths 
through the medium, the energy loss acquires a probabil- 
ity distribution which is then convoluted with the pQCD 
parton spectrum [25l . [2^. A direct comparison [77| of 
model results with a computation retaining the dynami- 
cal fluctuations suggests that the tomographical interpre- 
tation of the result is drastically different. While in an 
average picture parton propagation through the medium 
center is strongly suppressed, strong fluctuations make it 
a likely scenario. 

If one considers the strength of the fluctuaions, there 
is no compelling reason to assume that dynamical fluc- 
tuations can be averaged out. Plots of P{AE) in radia- 
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tive energy loss formalisms typically do now show a nar- 
row peak around the mean energy loss, but a wide, non- 
Gaussian distribution probing the whole available range 
of energies |78|. For incoherent processes, this has been 
also demonstrated in a realistic treatment of elastic en- 
ergy loss [H^l where it was also found that the dynamical 
fluctuations are strong and non-Gaussian. 



4- Implementation details 

Different implementations of the same physics idea 
can differ in parameters such as cutoffs, regulators to 
treat kinematical situations outside the model assump- 
tions and similar details. A very instructive review com- 
paring the ASW, GLV, AMY and HT formulations of 
radiative energy loss in a fixed-length constant de nsity 
medium ('brick') illustrates the effect of these details |78j . 
and wc refer the reader to this work for an in-depth dis- 
cussion of these effects. To summarize the findings very 
briefly: In an extreme case, the same observed suppres- 
sion pattern would lead to a factor 8.5 (!) difference in 
the extraction of the medium transport coefficient de- 
pendent on which formulation is used. This uncertainty 
is largely driven by implementation-specific details and 
hence is present in addition to all the uncertainties not 
related to the parton-medium interaction as discussed in 
the previous section. 



tions) indicates that the approximations may be funda- 
mentally fiawed and too simplistic. However, these un- 
certainties mainly seem to influence the overall strength 
of the parton-medium interaction. In (28| it is demon- 
strated for ASW, AMY and HT that despite the huge 
discrepancies in the extracted medium parameter, once 
the interaction strength is fixed for each model to the 
data in one place, the resulting dynamics in terms of Pt 
and centrality dependence or even unobservable quanti- 
ties like hadron escape probability are similar to the 10% 
level. 

This, taken together with the observations of the pre- 
ceding section that the uncertainties associated with a 
particular model assumption affect various observables in 
a very different way, gives rise to the hope that by test- 
ing models against a specific set of observables which are 
dominantly sensitive to a generic feature of models, one 
can constrain the physics of jet quenching in a meaningful 
way without overly being affected by specific implemen- 
tation details. In other words, an incoherent model of 
parton-medium interaction will show a linear pathlength 
dependence of energy loss under very general conditions, 
quite independent of implementation details. If it can be 
demonstrated that a particular observable probes path- 
length dependence and rules out linear dependence, then 
a whole class of models can be ruled out without refer- 
ence to specific details. 



5. Summary 



IV. CONSTRAINING PARTON-MEDIUM 
INTERACTION PHYSICS 



The sensitivity to assumptions discussed in this section 
fall into rather different categories. Some are directly 
linked to research questions, for instance the mass of the 
medium constituents probed by a hard parton which reg- 
ulates the relative strength of elastic energy loss is not 
just a model parameter, but a probe characterizing the 
mircoscopical degrees of freedom in the medium, i.e. a 
quantity to be determined by comparison with experi- 
ment. Qualitatively different dynamics based on e.g. the 
choice of the parton-medium interaction is best not dis- 
cussed in terms of uncertainties but in terms of different 
a priori plausible scenarios to be compared with data. 

The sensitivity to other assumptions should not be seen 
as uncertainty in the modelling — if an approximation 
of constant energy loss given a path results in a quali- 
tatively different picture of the process than a more de- 
tailed computation allowing for fluctuating energy loss 
around its mean value, then this demonstrates just that 
the approximation is not justified. Similarly, wherever 
a full in-medium shower evolution shows qualitative dif- 
ferences to an energy loss calculation, this signals the 
breakdown of the energy loss approximation rather than 
an uncertainty in the calculation. 

The sensitivity to implementation specific details as 
demonstrated in (TSj is most troublesome, as this (in the 
absence of a clear picture of a hierarchy in the assump- 



In the following, the aim is to discuss constraints from a 
set of observables, each of which is in particular sensitive 
to few or ideally one generic model property. Thus, each 
of the observables provides different qualitatively differ- 
ent constraints which can be used to rule out classes of 
model with certain properties. In the following we will 
argue that the data on Raa{Pt, </', a/s) Iaa{pt, \fs) 
constitutes such a set of observables. 

Among these observables, Raa{Pt) at RHIC kinemat- 
ics exhibits a fairly generic Pt dependence, but can be 
used as a starting point to fix all model parameters given 
the choice of a particular hydrodynamical model. The 
dependence of Raa on the reaction plane angle (j), in 
particular the spread S^^^ between emission into the re- 
action plane and out of the reaction plane is then, to 
almost equal parts, sensitive to the choice of the hydro- 
dynamical model and to the pathlength dependence of 
the parton-medium interaction. The requirement that 
the same model describes Raa{Pt) at LHC energies of 
2.76 ATcV when a carefully constrained extrapolation of 
the medium model to larger ^/s is done then probes the 
parton-medium interaction model-specific probability to 
find an almost unmodified fragmentation pattern. Fi- 
nally, Iaa can be used to observe the more differential 
pattern of energy redistribution for subleading shower 
hadrons once the pathlength dependence is fixed. 
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It needs to be stressed again that the large medium 
evolution model dependence of observables makes it 
mandatory to use a medium evolution model which is 
constrained by bulk data. 



A. Pt dependence of Raa at RHIC 

The Pt dependence of 7r° Raa in central 200 AGeV 
AuAu collisions as obtained by the PHENIX collabora- 
tion [t^ is usually the first high Pt observable to be 
discussed in models. 

For a wide range of model assumptions, the Pt de- 
pendence of Raa in this kinematical situation is generic: 
In [3^ it is shown that no large difference between the 
ASW, HT and AMY formulation of energy loss can be 
seen and in [2^, Is^l no visible dependence on the assump- 
tions about medium modelling are observed. An expla- 
nation for this finding is provided in (76j by an analysis 
how sensitive the convolution of parton spectrum and en- 
ergy loss probability density P{AE) is to the choice of 
P{AE). 

The main findings can be illustrated in a simple model. 
Consider suppression in the leading parton energy loss 
picture on the partonic level only. P{AE) can then have 
a discrete transmission piece T representing the proba- 
bility for no energy loss and a continuous piece P'{AE) 
representing a parton shift in energy by AE. After the 
convolution with a steeply falling parton spectrum, this 
is equivalent to a suppression since a parton which is after 
the energy shift found at pt must have originally come 
from a higher momentum pT + AE where less partons 
are available. Here it is understood that AE > E implies 
that the parton is absorbed into the thermal medium. 

Raa can then be understood from the ratio of modi- 
fied over unmodified parton spectrum, where the modi- 
fied parton spectrum at a given pT is determined by the 
number of partons escaping without energy loss plus the 
number of partons available in the spectrum at pt + AE 
times the probability P{AE) of a shift by AE. If we 
assume a power law p^" for the parton spectrum, 



Raa{pt) ~ T ' 



d{AE)P'{AE) 1 



AE\ 
PT J 



(10) 

It is evident from the expression that Raa at a given pT 
is equal to the transmission term of zero energy loss plus 
a contribution which is proportional to the integral of 
{P{AE))taa from zero up to the kinematic limit energy 
scale E'rnax -56671 through the filter of the steeply falling 
spectrum. Raa then generically grows with pT since 



the 'penalty factor' 



1 + 



PT I 



for observing a parton 



downshifted by AE from a higher momentum scale de- 
creases with increasing pT (in reality, the last statement 
is no longer true when the actual pQCD spectrum de- 
viates significantly from a power law, i.e. around -s/s/^i 
however for realistic experimental conditions this region 
cannot be accessed with meaningful statistics). 



The speed of growth depends on the weight of 
{P{AE))taa ™- the region from zero to iSmax and on the 
power n of the parton spectrum. 

At RHIC kinematics, the spectral power is n ^ 7.. 8 
and in the limited kinematic range between 6 GeV where 
a perturbative description of the spectrum becomes avail- 
able and 15 GeV where the experimental errors become 
substantial, a (moderate) 5 GeV shift in parton energy al- 
ready corresponds to 80% suppression. This means that 
the Pt dependence of Raa in this domain only probes 
the probability that an energy loss model results in shifts 
between zero and 3.. 4 GeV. In contrast, energy loss 
models typically expect shift probabilities throughout the 
whole allowed kinematic range O(IOO) GeV [zi,!^]- The 
implication is that the experiment is simply not sensitive 
to a large region of P{AE), instead fixing the normal- 
ization of models to data already limits the possible Pt 
dependence to a range which is not distinguishable from 
a constant value given the experimental errors. 

This means that any model which can be described ef- 
fectively in terms of P{AE) result in approximately the 
same Pt dependence of Raa, no matter what model as- 
sumptions are made in detail. The obvious exception are 
models which have an explicit dependence on the initial 
parton energy E in the parton-medium interaction, i.e. 
can be effectively cast into the form P{AE\E). At least 
two such scenarios have been investigated: In [76j it was 
established that a fractional energy loss AE = fE leads 
to a decrease of Raa with Pt which is not supported by 
the data. The resummed HT formalism [s^ and YaJEM- 
D (6l| also contain, as discussed above, such an explicit 
dependence which in this case leads to a stronger growth 
of Raa with Pt than in other models. 

It follows from the above that the only really non- 
trivial information from Raa{Pt) in central 200 AGeV 
collisions is the absolute normalization, i.e. the overall 
strength of the parton-medium interaction. Thus, given 
a chosen combination of medium evolution and parton 
medium interaction model, this data set can be used to 
e.g. fix the relation between thermodynamical parame- 
ters and transport coefficients, or to determine any other 
model-specific parameter which regulates the interaction 
strength. Once this is done, the combination of models 
to be tested is typically fixed for ^/s ~ 200 AGeV, i.e. 
different centralities at the same collision energy or dif- 
ferent colliding ions can be computed without additional 
free parameters. 



Reaction plane angle dependence of Raa at 
RHIC 



The first non-trivial test of a parton-medium interac- 
tion model is then if it predicts the centrality depen- 
dence of Raa correctly. Note that hydrodynamical mod- 
els are typically constrained using V2 and bulk Pt spec- 
tra across various centralities, thus while medium density 
and geometry changes from central to peripheral colli- 
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sions, these changes happen within 'the same' hydrody- 
namical framework and are not arbitrary. In particular, 
they do not usually involve additional parameters which 
would affect the high Pr side of the modelling. 

Initially it was assumed that predicting the central- 
ity dependence of Raa{Pt) averaged over reaction plane 
angle would be a reasonable test for models, and con- 
sequently this has been discussed in e.g. [s^. However, 
note that the two endpoints are already constrained — in 
central collisions, Raa is described by the models because 
this is the point where model parameters are adjusted 
(see above), Raa in ultra-peripheral collisions needs to 
go to unity because such collisions can not be different 
from a p-p collision, and models must monotonically con- 
nect these two limits, since both parameters relevant for 
the strength of energy loss, i.e. mean in-medium path- 
length and mean density arc monotonically decreasing 
from central to peripheral collisions. This means that 
studying angular averaged Raa is not a particularly sen- 
sitive test of models. 

This is different if the angular dependence of Raa is 
taken into account and the suppression is studied as a 
function of the angle of a high Pt hadron with the 
bulk matter reaction plane. For events inside a given 
centrality class, neither the mean density of the bulk 
medium nor the mean geometry changes, but the mean 
pathlength which a particle spends in matter is very dif- 
ferent between in-plane (i.e. along the short side of the 
almond-shaped overlap region) or out of plane (along the 
long side). Thus, the spread S'™^ betwen in plane and out 
of plane Raa is a direct probe of how parton-medium in- 
teraction models respond to a change in in-medium path- 
length. 

Of course, Sl^f docs not only depend on the pathlength 
dependence of the parton-medium interaction, but also 
crucially on how large the difference in mean pathlength 
between in-plane and out of plane emission is, and this 
in turn is often linked to the diffuscncss of the matter 
density profile — for instance, in [81j], striking differences 
between a hard sphere and Gaussian overlap models were 
found. 

Embedding a high pT event into a hydrodynamical 
bulk medium description never probes pathlength L with 
the same power as in the constant medium case. The 
longitudinal Bjorken (e.g. [l^) or close to Bjorken (e.g. 
(20j ) expansion implies that medium density drops like 
~ 1/r (where r is the proper medium evolution time) 
and this effectively removes one power from the path- 
length dependence. The transverse expansion velocity 
of the medium is initially small and grows over time as 
transverse pressure gradients lead to the development of 
flow. The resulting transverse expansion causes two po- 
tentially opposing effects — as the medium expands, the 
distance to the relevant decoupling surface, and hence 
the in-medium pathlength may grow as compared to a 
trans versally non-expanding case, whereas the medium 
density is decreased. In the case of a trans versally ex- 
panding cylinder and dependence of energy loss, these 



two effects would exacly cancel. However, in any real sys- 
tem the dynamics depends on details — the decoupling 
surface often moves inward in spite of transverse expan- 
sion, the rate at which density is diluted transversally 
depends on the stiffness of the equation of state, i.e. how 
efhciently spatial gradients are mapped into momentum 
space, and fluctuations in the initial state introduce fur- 
ther small-scale gradients and complicate the problem. 
Thus, there is no simple analytic way to understand how 
the mean pathlength difference between in-plane and out 
of plane is generated in a given hydrodynamical evolu- 
tion. 

For a representative set of four different hydrodynam- 
ical models, this problem has been studied in [s^l- Here, 
it was found that what matters most is late time dynam- 
ics in terms of how far out in terms of space and time the 
medium geometry extends, followed by initial dynamics, 
viscosity and initial profile diffuseness. The following ta- 



summarizes the findings: 


property of hydro 


rel. sensitivity So"t 


spatial size of freeze-out surface 


~ 50% 


initial time scale 


~ 25% 


viscous entropy generation 


~ 18% 


initial profile diffuseness 


~ 7% 



Thus, contrary to intuition, the actual choice of the 
initial matter density profile diffuseness is not a substan- 
tial effect. However, given this large variation of mean 
in-plane and out of plane pathlength between medium 
models which translates into a sensitivity of 5'™j to the 
choice of a medium model which is as large as the sensi- 
tivity to the choice of a parton-medium interaction model 
[soj , the implication is that pathlength dependence of en- 
ergy loss can not be studied in a meaningful way unless 
the full systematics on the hydrodynamical side is ac- 
counted for. 

As indicated above, the change of angular-averaged 
Raa as a function of centrality is constrained at the 
two endpoints of central and ultra-peripheral collisions. 
However, the evolution between these centralities de- 
pends both on the change in mean density and mean 
pahtlength, thus one would expect an ordering of models 
dependent on what power of the pathlength is probed 
by models, with the incohcrend, L-depcndcnt models re- 
sponding slowest to a change in centrality and the AdS- 
type strong coupling dependent models responding 
relatively faster. Such trends are indeed observed when 
compa ring par ton-medium interaction models which each 
other [30, [S^l however except in the case of linear 
pahtlength dependence they are not very significant with 
the data. 

The combined results of [13, HH, HI, H!] , i.e. the viabil- 
ity of various physics mechanisms underlying a particular 
pathlength dependence given a choice of medium model 
and the data is best summarized in the form of a matrix 
(Tab.m. 
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model 


elastic L radiative 


AdS 


rad. fin. 


E min. Qo 


3+ld ideal 


fails 


works 


fails 


fails 


works 


2+ld ideal 


fails 


fails 


marginal 


fails 


fails 


2+ld vCGC 


fails 


marginal 


works 


fails 


marginal 


2+ld vGlb 


fails 


marginal 


works 


fails 


marginal 



TABLE I: Viability of various combinations of medium model and parton-medium interaction model given the PHENIX data 
for Raa{4>) [z3 (based on calculations in [sol. [6ll. Is^. Is^ . see text for the meaning of the various labels for models). 



It can be observed that two classes of models fail un- 
conditionally, i.e. regardless what hydrodynamical model 
is chosen, i.e. the elastic (or incoherent) L-dependent 
type and the radiative type with finite energy corrections 
which for longer path is effectively L-dependent. The fail- 
ure is quite dramatic — S'™^ differs from the data by at 
least a factor 6. For all other model types, it is possible 
to find a hydrodynamical evolution which results in a fair 
description of the data. 



AuAu 200 AGeV, 40 - 50 % 
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FIG. 2: Comparison of in-plane and out-of plane Raa for an 
elastic MC model [sl, YaEJM [6l[ and the ASW radiative 
energy loss model [5J] using the same 2-1- Id hydrodynami- 
cal medium evolution [31 as compared with 40-50% central 
200AGeV AuAu collision data (circles: in plane, diamonds: 
out of plane) obtained by PHENIX 



given the data. This is perhaps not surprising, as it has 
never been suggested that coherence is completely ab- 
sent. However, given the fact that no model overpredicts 
^out by more than 15%, it can be inferred that an inco- 
herent component must be ^ 10% in order for a combined 
incoherent -I- other mechanism scenario to describe the 
data [H,!!!]. 

This is in contrast to straightforward computations of 
the relative magnitude of elastic energy loss from pQCD 
where f req uently elastic components 0(50%) are found 
13 US HJ] ■ The inevitable conclusion is that the main 
assumption made in pQCD calculations of elastic energy 
loss, i.e. that the medium DOF are almost free (quasi)- 
particles which can take a sizeable amount of recoil en- 
ergy away from a leading parton does not appear to be 
true in nature. 

It is also very instructive to note that YaJEM in which 
the linear pathlength dependence arises only effectively 
due to finite energy and length corrections to a medium- 
induced radiation formalism which takes coherence into 
account and which in principle exhibits dependence 
fails just in the same way as a formalism which is inco- 
herent from the beginning [6l| . This indicates that the 
observable really probes pathlength dependence regard- 
less of how it arises microscopically. The fact that finite 
energy corrections lead to the failure of a radiative energy 
loss formulation seriously questions the relative success 
of radiative energy loss models derived in infinite kine- 
matics (such as ASW or the AdS model). However, this 
is a qualitative theoretical objection and as we shall see 
in the next section, these models can also be ruled out 
based on a comparison with data. 



These findings are illustrated in Fig. [2] at the example 
of the elastic energy loss MC calculation [H], YaJEM 
(el [72! and the ASW radiative energy loss model ^ 
embedded into a 2-|-Id ideal hydrodynamical model for 
the bulk medium evolution. With this particular medium 
model, radiative energy loss misses S'™^ by about a factor 
two, but clearly reaches the right order of magnitude and 
hence is compatible with the data if the full uncertainty 
associated with the medium model choice is taken into 
account. In contrast, both the elastic model and YaJEM 
fail clearly in reproducing the spread and the average 
even allowing for the systematic uncertainty associated 
with the choice of the medium evolution. 

The immediate conclusion is that a purely incoher- 
ent parton-medium interaction picture can be ruled out 



C. Pt dependence of Raa at the LHC 

Following the arguments given in IIV A[ we can con- 
clude that at LHC kinematics the harder primary parton 
spectrum (i.e. the lower power n = 4..5 in a power- 
law fit) will allow to observe a larger range of P{AE) 
and that thus the Pt dependence of Raa in central 2.76 
AGcV central PbPb collisions will probe, unlike in the 
RHIC case, model-specific differences. 

If one aims at a comparison of a model combination 
which is constrained by RHIC data with Raa at LHC, 
there are however a number of additional difficulties con- 
nected with extrapolating a hydrodynamical framework 
to higher -y/s (for an extended discussion see [l^). For 
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instance, it is expected that for increasing -y/s the parti- 
cle production mechanism might change soft participant 
scaling for the produced entropy at low ^/s to hard bi- 
nary collision scaling at higher ^/s, thus even qualita- 
tively changing the initial matter profile. Likewise, the 
matter equilibration time is expected to shorten, while at 
the same time the overall higher multiplicity would allow 
for a lower decoupling temperature for the bulk matter. 
Thus, unlike in the previous extrapolation to different 
centralities in the same colliding system, there is some 
reason to expect that additional free parameters may ap- 
pear in the extrapolation of the hydrodynamical medium 
model from RHIC to LHC energies, some of which may 
be relevant for high Pt physics as well. It is therefore 
unreasonable to expect that a RHIC constrained model 
would exactly match LHC data if the 'same' hydrody- 
namical model is used at both RHIC and LHC, since 
the definition of 'same' on the hydrodynamical side now 
hinges on details of the framework chosen to extrapolate 
the hydrodynamics. In the following, we employ a pQCD 
-I- saturation approach, the "EKRT model" [8g as the 
relevant framework to define the 'same' hydrodynamics 
for different ^/s fssj . 

Note that this specifically requires the use of the 
2-1- Id hydrodynamical framework [l^ as the underlying 
medium model which works worst with many parton- 
medium interaction models for angular differential ob- 
servables at RHIC energies (see llV Bl above). Thus, while 
we know that Raa{Pt) for central collisions is not very 
sensitive to the geometry of the hydrodynamical mod- 
elling, there is no reason to expect that this calculation 
agrees well with angular differential data (since it mani- 
festly does not at RHIC). 

Thus, what is tested by this investigation is the re- 
sponse of parton-medium interaction model to a signifi- 
cant change in medium temperature (or equivalently den- 
sity of scattering centers) and to the hardening of the 
primary parton production pT spectrum. While the first 
effect primarily governs an overall normalization factor 
of Raa since the strength of parton-medium interaction 
scales with the density of scatterers, the second effect 
leads to explicit Pt dependence. 

A comparison of various parton-medium interaction 
model when extrapolated from RHIC to LHC energies 
using the closed EKRT framework is shown in Fig. [3] (see 
[ssj l for details of the calculation). It is evident that the 
Pt dependence of various models indeed differs as ex- 
pected, and also that the normalization of a suitable Pt 
average is quite different, i.e. we observe both the re- 
sponse to the changed medium density and to the hard- 
ening of the parton spectrum. 

It is clearly seen that both the ASW and the AdS 
frameworks (and to a lesser degree YaJEM) 'ovcrqucnch', 
i.e. they predict more suppression than seen in the data. 
The overquenching of the ASW model appears to be 
generic for radiative energy loss models, as it also appears 
for YaJEM and a qualitatively rather similar result was 
obtained within the WHDG model 
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FIG. 3: Various RHIC-constrained parton-medium interac- 
tion models (see text) extrapolated to 2.76 ATeV central 
Pb-Pb collisions at LHC using the EKRT minijet saturation 
model [s^ to extrapolate a 2-|-ld hydrodynamical evolution 
along with a pQCD calculation for hard parton production 
[sst compared with ALICE data 



Given the overquenching of radiative energy loss mod- 
els, an even stronger overestimation of the parton- 
medium interaction strength in a strong coupling descrip- 
tion (AdS) is inevitable: As seen in IIIIB 11 while para- 
metrically the energy loss in radiative models scales as 
^ T'^L^, in a strong coupling description this becomes 
^ 2"4^3 gjjjcg a,t higher ^/s both mean temperature and 
(to a lesser degree) mean size of the medium increase, 
paramctrically Raa in a strong coupling description must 
extrapolate in ^ to a lower value than in a radiative de- 
scription, thus if radiative models overquench, a strong 
coupling description will fare even worse. 

With regard to the Pt dependence, one can observe 
that it appears to be too weak in the majority of models, 
implying that the probability to loose a small amount of 
energy from the leading parton is underestimated by such 
models. The exception to this rules are parametrized 
elastic models with a Gaussian assumption for P{AE) 
[sst which by construction leave more room for small en- 
ergy losses as compared with the almost fiat P{AE) char- 
acteristic for radiative energy loss models, and YaJEM-D 
which has an explicit energy dependence, i.e. can only 
be cast into the form P{AE\E). Since we already ex- 
cluded a dominant contribution of an incoherent energy 
loss mechanism based on pathlcngth dependence in the 
previous section, the only relevant model at this point 
is YaJEM-D. In comparison with YaJEM, it can be ob- 
served that this explicit reduction of energy loss with in- 
creasing initial parton energy E leads both to a different 
normalization and a steeper increase with Pt, in much 
better agreement with the data. 

The same mechanism as in YaJEM-D (i.e. minimum 
scale down to which a shower is evolved in the medium 
determined by Qq = \/ E / L) has been used, albeit in a 
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simplified implementation, in the context of the HT for- 
malism [s^ , and qualitatively in agreement with YaJEM- 
D, a strong increase of Raa with Pt and no ovcrquench- 
ing is observed. 

We can infer from these observations that Pt, ^/s de- 
pendence of Raa and centrality dependence provide in 
some sense orthogonal constraints. While the centrality 
dependence favours models in which energy loss scales 
with a high power of pathlength and rules out incoherent 
mechanisms, almost the opposite is true for the Pt de- 
pendence at LHC where models with a weak pathlength 
dependence fare better with the data since they have a 
higher probability of inducing only small energy loss. 

The only class of models which can account for both 
observables simultameously introduces a minimum in- 
medium virtuality scale, as exemplified by both YaJEM- 
D [6l| and resummed HT [5^ [s^. This implies that 
the dominant mechanism of leading parton energy loss is 
indeed likely to be medium-induced perturbatively cal- 
culable radiation, but this mechanism can not be treated 
in a leading parton energy loss model, since such a model 
can not implement a constraint on the minimum virtu- 
ality scale of a developing shower. Thus, coherence is 
an important ingredient, but on the level of the LPM 
effect as implemented in leading parton energy loss mod- 
els it is effectively not much different from an incoherent 
mechanism j6(]| - l62| . Currently, there is no evidence that 
a strong coupling scenario is favoured (or even allowed) 
by the combined high Pt observables. 



D. zt dependence of Iaa at RHIC 

The suppression factor Iaa of back-to-back correlated 
yield corresponds to a conditional probability of finding 
subleading hadrons on the near or away side given a trig- 
ger hadron in a certain momentum range. The require- 
ment of a trigger biases the selection of events in various 
interesting ways as compared with a full jet observable 
where (in principle) every hi gh p t scattering process con- 
tributes in an unbiased way [90[ . With zt being the ratio 
of trigger momentum over associate hadron momentum, 
the low Zt range of Iaa then probes the fragmentation 
pattern of subleading shower hadrons given the biased 
selection of triggered events. 

Due to the combination of biases, Iaa is a fairly com- 
plicated observable, but it has several advantages over 
single inclusive hadron suppression. First, since the path 
of the away side parton through the medium is always 
different from the path of the near side parton, a simul- 
taneous measurement of Raa and Iaa allows to probe 
pathlength dependence of a parton-medium interaction 
model for central collisions, i.e. without any change in the 
hydrodynamical evolution from the point where model 
parameters are calibrated. Second, the near and away 
side correlated yield at low zt is sensitive to subleading 
fragmentation, i.e. probes the dynamics of energy redis- 
tribution within a shower beyond the validity of a leading 



parton energy loss model. Third, the trigger requirement 
leads (for RHIC kinematics more strongly than for LHC 
kinematics) to a bias for the away side parton to be a 
gluon and thus opens some sensitivity to the parton type. 

The systematics of the dependence of Iaa both on the 
medium evolution model and the parton-medium inter- 
action model has been investigated in [oOl and [U . 

The main results can be summarized as follows: For 
Zt > 0.5 (i.e. for leading hadron physics), the sys- 
tematics of Iaa agrees with what can be obtained from 
Raa{4')-: the same pathlength dependences in com- 
bination with the same hydrodynamical models are pre- 
ferred by both observables, although there are differences 
in the sensitivity in detail. 

For Zt < 0.5, clear differences between models uti- 
lizing the energy loss approximation and full in-medium 
shower evolution codes are apparent. While energy loss 
models qualitatively show the wrong behaviour with the 
data, i.e. a decrease of Iaa when approaching low zt 
from above, in-medium shower codes show an increase 
corresponding to the medium-induced additional gluon 
radiation. Note that in energy loss models this contri- 
bution is absent because the fate of subleading partons 
is not explicitly tracked by definition, and thus energy 
lost from the leading partons by construction disappears 
from all observables. 

AuAu 200 AGeV 0-5% centrality 

triggers - 15 GeV 
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• YaJEM-DE 
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FIG. 4: Away side dihadron suppression factor Iaa computed 
with YaJEM in a pure radiative scenario (YaJEM-D) and 
with a 10% contribution of elastic energy loss (YaJEM-DE) 
[qH and for comparison also in the energy loss approximations 
models ASW and AdS using a 3-1-1 d hydrodynamical model 
for the medium [S^l compared with STAR data [o^l . Lines in 
theory results are drawn to guide the eye. 

The observation of an upward trend of Iaa for low 
Zt in itself can thus be taken as evidence for a medium- 
induced radiation scenario, as this region probes the rem- 
nant of the radiation after hadronization explicitly. This 
is illustrated in Fig. U with results from YaJEM-D us- 
ing the best choice for the medium model in comparison 
with model results using the energy loss approximation. 
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Clearly, the two energy loss models ASW and AdS miss 
the trend visible in the low zj- data. As a side remark, 
note that the AdS model as shown in the figure overpre- 
dicts the amount of suppression seen in the data even 
at high zt- This is driven by the choice of the hydro- 
dynamical model, which is here selected to give the best 
fit using YaJEM-D. Based on such a figure alone, one 
might be tempted to rule out AdS as not viable with the 
data, however if a different medium model ([3 i^" 
stance) would have been chosen, AdS would agree well 
with this set of data whereas YaJEM-D would not. This 
strongly emphasizes the point that the systematics asso- 
ciated with the medium modelling must be understood 
properly before drawing conclusions. 

One can observe that a fully radiative scenario as ex- 
emplified by YaJEM-D overpredicts the rise at low zt- 
However, any elastic contribution to the parton-medium 
interaction can be expected to deplete the medium- 
induced radiation spectrum by transfering energy from 
soft gluons into the medium as suggested e.g. in [93| . 
Allowing for a 10% elastic contribution (i.e. about the 
maximal value that is allowed by pathlength-dependent 
observables) good agreement with the data is obtained 
[qH . Note that such a contribution is small enough not 
to change any of the single hadron observables discussed 
above in a significant way, i.e. such a scenario (referred 
to as YaJEM-DE in the following) agrees with the results 
shown for YaJEM-D above. 

We can infer from these results that the zt dependence 
of Iaa confirms the constraints for pathlength and geom- 
etry as outlined above. In addition, the observation of the 
onset of the medium-induced radiation spectrum argues 
for a dominantly radiative energy loss picture and makes 
the treatment in terms of in-medium shower evolution 
mandatory — energy loss modelling leads to even quali- 
tative failure to describe the data at low zt- The amount 
of medium-induced radiation is overpredicted by a pure 
radiative picture and hence allows to constrain the con- 
tribution of elastic interactions from below. This comple- 
ments the constraints on elastic interactions from above 
by pathlength-sensitive observables, and both sets of con- 
straints can be satisfied with a fraction of ^ 10% elastic 
energy transfer. In principle, this value contains valuable 
information about the microscopic dynamics of parton- 
medium scattering (for instance, the amount of elastic 
energy transfer is a function of quasiparticle mass). 



E. Summary 

In summary, a combined analysis of models in the con- 
text of the measured Raa{Pt, 4>, \/s and lAAizr), taking 
into account the systematic uncertainties associated with 
the medium evolution model suggests the following pic- 
ture: 

• The dominant fraction of the parton-medium in- 
teraction is coherent. A large incoherent contri- 



bution to the interaction necessarily leads to a lin- 
ear parthlength dependence of energy loss, and this 
can be ruled out with great confidence due to the 
failure of linear pathlength models to describe the 
reaction plane dependence of Raa- These findings 
are confirmed by the failure of linear pathlength 
dependence to account for Iaa [83| . 

• Pathlength dependent observables rule out in addi- 
tion to incoherent mechanisms also models in which 
coherence playes a role, but in which corrections 
such as finite kinematics or finite length effects ef- 
fectively reduce the pathlength dependence to lin- 
ear. Several theoretical results |60l - [63 | suggest that 
this happens generically for radiative models im- 
plementing the LPM effect which generate an 
dependence of energy loss for infinite kinematics, 
i.e. the constraints from pathlength dependence 
are quite severe and affect a large class of suggested 
scenarios. 

• While S'™j tends to favour the scaling of energy loss 
with a high power of pathlength, dimensional argu- 
ments then suggest that energy loss then also scales 
with correspondingly higher power of the medium 
temperature T. Such a scaling is not indicated by 
the observed Raa at the LHC — radiative energy 
loss models parametrically scaling with tend to 
overpredict the amount of suppression, correspond- 
ingly any model scaling energy loss ~ L" with 
n > 2 will scale even worse as long as the average 
temperature in the LHC medium is higher than in 
the RHIC medium. This disfavours strong coupling 
models which can generate scaling with higher pow- 
ers of L. 

• The combinaion of the above two points suggests 
that the pathlength dependence of energy loss is 
not just a simple power of T as given e.g. by the 
LPM effect assuming infinite parent parton kine- 
matics. This is also suggested by more detailed the- 
oretical arguments. One scenario that is allowed by 
the data is a restriction on the minimum virtuality 
in the shower that can be probed in the medium as 
used by YaJEM-D or the resummed HT approach, 
other scenarios may be possible. 

• The observation of a rise in the low zt region of 
Iaa is consistent with a radiative energy loss sce- 
nario. Qualitatively, additional medium-induced 
radiation is expected to produce such a rise. Quan- 
titatively, a constrained calculation with a small 
fraction of elastic energy loss (YaJEM-DE [9l|) is 
able to account for the data. In this framework, the 
elastic fraction of energy loss can be constrained 
from above and from below to be about 10%. The 
data at low zt cannot even qualitatively be ex- 
plained in an energy loss picture, a full treatment 
of subleading shower partons is needed. 
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RAA{<l>)@RmC Raa@LRC (Pt) /aa@RHIC 


elastic 


fails 


works 


fails 


elMC 


fails 


fails 


fails 


ASW 


works 


fails 


marginal 


AdS 


works 


fails 


marginal 


YaJEM 


fails 


fails 


fails 


YaJEM-D 


works 


works 


marginal 


YaJEM-DE 


works 


works 


works 



TABLE II: Viability of different parton-medium interaction 
models tuned to the Pt dependence of Raa in 200 AGeV Au- 
Au collisions given various data sets under the assumption 
that the best possible hydrodynamical evolution scenario is 
chosen. The various labels refer to: elastic fs^. elMC [4^ . 
ASWjp, AdS [H, YaJEM [zl.lzl, YaJEM-D^, YaJEM- 
DE 



If we always select the optimal medium evolution 
model for a given parton-medium interaction model, the 
findings of the systematic analysis can be summarized in 
table ini 

It becomes readily apparent that some constrains work 
in an orthogonal way, for instance while Raa at LHC 
does not lead to conclusive statements about elastic en- 
ergy loss scenarios, any pathlength dependent observable 
quickly establishes their failure. 

Of the models for which the full systematics has been 
tested, only one combination (YaJEM-DE ^l| in the 
34-ld ideal hydrodynamical code by Nonaka and Bass 
[20j results in a viable description of all the combined 
data discussed so far. Other models (for instance the 
HT approach) for which the full systematics with realis- 
tic medium models is not known exhibit the properties 
required by the data, and there is no reason to assume 
such models would not lead to an equally good descrip- 
tion of the data. Given that the full matrix of seven a 
priori viable parton-medium interaction models (based 
on a description of Raa{Pt) in 200 AGeV central Au-Au 
collisions at RHIC) combined with four possible choices 
of the medium evolution model for which the systemat- 
ics is known contains 28 different combinations, the con- 
straining power of the multi-observable analysis becomes 
readily apparent. 



V. OPEN ISSUES 

This section summarizes various observables or model 
frameworks which have, for reasons given for each in the 
specific context, not been included into the analysis of 
the previous section. 



A. Raa of protons 



As discussed in IIII B 21 one of the central assump- 
tions underlying the theory of hard probes is that the 



final state interaction is predominantly a partonic phe- 
nomenon, i.e. that the observed effect is not driven by 
hadron-hadron interactions during the later stage of the 
medium evolution. The obvious way to test this assump- 
tion is to study observables for various identified hadron 
species. 

If the final state interaction is predominantly a par- 
tonic phenomenon, then Raa for various hadron species 
should depend on the relevant partonic production chan- 
nels for the hadron, but not on the actual identity of the 
hadron. Within the AKK set of fragmentation functions 
jo^ the proton is found to be dominated by the frag- 
mentation of gluons in the RHIC kinematic range 5-10 
GeV whereas the pion receives about equal contribution 
from quark and gluon fragmentation. Since gluons in- 
teract more strongly with the medium by a color fac- 
tor Cp — 9/4, the natural expectation is that proton 
R^j^{Pt) < R^aIPt)- Experimentally however, the op- 
posite trend was observed [9a |. 

On face value, this finding argues against a partonic 
origin of the suppression and it is currently not well ex- 
plained. In a detailed computation using the ASW for- 
malism, the nuclear suppression factor of pious and pro- 
tons is found to be fairly similar, but as expected with 
R\j^{Pt) < R\j^{Pt) |9g|. If conversion reactions like 
qq — )■ gg which change the identity of the hard parton 
are introduced [o^l then proton and pion Raa can at 
best be made equal. It should be noted however that one 
of the central assumptions made in almost every model 
of parton-medium interaction, i.e. that the hadroniza- 
tion process takes place outside the medium, is violated 
in the kinematic range of the data. If the formation time 
of protons is estimated as r ~ Ep/rUp a value of ^ 2 fm is 
obtained, i.e. it can not be argued that proton formation 
happens outside the medium, instead the scale indicates 
that in most cases it happens inside the medium. Thus, 
models which assume a vacuum hadronization pattern 
are not applicable to this data. 



B. Heavy quark energy loss 

The idea to study the energy loss of heavy c and b 
quarks is based on the notion that their radiative energy 
loss contribution should be significantly reduced as com- 
pared to light quarks due to the so-called dead cone effect 
[98| . This has led to the expectation that heavy quark 
Raa (as probed by single inclusive electron suppression) 
should be observed above light quark driven hadronic 
Raa- However, the single inclusive non-photonic electron 
suppression measurements (which predominantly probes 
heavy quark decays) have shown a comparable magni- 
tude of heavy and light quark suppression [o^ Il00| . This 
has sometimes been refered to as 'heavy quark puzzle' 
and lead to the coiyecture of a large elastic energy loss 
component ^ 50% \5l\ . 

Such a large component in the light quark sector is 
however, as we have seen in previous sections, con- 
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clusively ruled out by pathlength-sensitive observables. 
Thus, one is forced to either conclude that a large elastic 
component is not the solution to the heavy quark puzzle 
or that the physics of light and heavy quark energy loss 
is significantly different. There is in fact some reason to 
suspect the latter scenario. 

First, the energy loss due to elastic collisions is ex- 
pected to be en hanc ed relative to light quarks due to the 
'tagging' effect |lOl| — while for an unidentified primary 
quark a collision leading to a transfer of 80% of its en- 
ergy to a medium parton effectively counts as 20% energy 
loss since the medium parton then becomes the new hard 
parton, a heavy quark is always identified since due to 
the high mass threshold it cannot be generated by ther- 
mal excitation and hence 80% energy transfer in this case 
effectively imply 80% energy loss. 

A second important aspect which is not obvious in 
an energy loss approximation becomes apparent in a 
medium-modified shower picture. The formation time 
of a heavy-quark induced shower, i.e. the time it takes 
to reduce the original quark virtuality at the hard scale 
to the quark mass scale, an d al so the formation time of 
heavy mesons is very short [l02j | . This would imply that 
the object interacting with the medium is i n fac t not 
a quark, but rather a pre-hadronic resonance jl02l Il03j 
which is quickly dissociated by the medium, but never- 
theless enhances the effective interaction cross section. 

Models based on elastic heavy- quark me dium interac- 
tion can explain the existing data (102| - |105 |. An alterna- 
tive approach is to postulate the same physics for light 
and heavy quark energy loss and treat heavy quark en- 
ergy loss using the strong coupling assumption [106[. As 
we have argued above, this is not in agreement with the 
scaling hadronic Raa from RHIC to LHC e nergi es in the 
light quark sector. As pointed out e.g. in jl07l |. expec- 
tations for heavy quark Raa at LHC are substantially 
different dependent on if a perturbativc or strong cou- 
pling scenario is used. 

In view of the fact that both the tagging and the for- 
mation time argument imply that the physics of heavy 
quark energy loss is different beyond a kinematic dead 
cone radiation suppression and may quite possibly con- 
tain non-perturbative QCD dynamics, the non-photonic 
single electron data has not been used for the analysis 
presented in this work. 



C. 7-hadron correlations 

In principle, 7-h correlations allow to constrain the 
energy of the recoiling parton at the hard vertex much 
better than dihadron correlations. In practice, the re- 
duced experimental statistics often requires a compara- 
tively wide trigger momentum window which tends to 
offset this advantage somewhat. While 7-h correlations 
probe the fragmentation function in the limit where the 
photon energy is known with great accuracy, the opposite 
situation in which photons within a wide energy range are 



accepted as triggers resembles more an Raa of quarks 
since the gq — >■ 7(7 production channel dominates over 
qq — > 79 and since, unlike for dihadron correlations, the 
photon does not bias the distribution of initial vertices 
given a trigger. 

From these considerations, it can be deduced that any 
parton-medium interaction model which is tuned to de- 
scribe single hadron suppression for a given choice of a 
medium mode l should a lso describe the high zt region of 
the 7-h data |l08l . Il09l | in the same medium model be- 
cause the integrals over geometry and energy loss prob- 
ability involved are just the same. This is indeed found 
for the Zhang-Owens- Wang- Wang (ZOWW) model [iToj . 
for the AMY model 11 1 ill for the ASW model and for the 
shower code YaJEM jll2l |. 

The interesting differences between the models (and 
between model and data) appear at low zt- For instance, 
as in the case of dihadron Iaa, YaJEM is found to over- 
shoot 7-h Iaa at low zt whereas ASW docs not show any 
rise coming from high zt down to low zt- This suggests 
the same basic physics interpretation: The rise at low 
Zt measures the strength with which medium-induced 
radiation leads to subleading shower hadrons. A com- 
putation of this contribution requires a full in-medium 
shower framework, and would constrain the contribution 
of elastic energy loss from below. 

At present, it is unclear if there is information in 7-h 
correlation that can not be obtained from the combina- 
tion of hadronic Raa and dihadron Iaa^zt)- For this 
reason and the lack of published systematic studies with 
different medium models, the observable has not been 
used for the analysis presented in this work. However, 
unlike in the case of proton Raa or heavy quark energy 
loss, there is currently no reason to assume that describ- 
ing 7-h correlations would be problematic for any model 
which fits the known constraints. 



D. The dijet asymmetry 

The energy imbalance between back-to-back jets in a 
medium has been quantified by ATLAS and CMS [11,111 
and has immediately sparked a lot of theoretical activity. 
Models able to account for the measured imbalance range 
from parametric estimates [93| via schematic modelling 
of energy loss and jet finding |113j to sophisticated NLO 
pQCD modelling of jet evolution in a medium, however 
without detailed medium or jet-finding modelling |114| 
and complete MC modelling of a jet embedded into a bulk 
medi um w ith realistic jet finding using the MARTINI 
code fTTsj . 

At this point, this clearly argues that there is noth- 
ing mysterious about the measured dijet imbalance which 
would require an exotic scenario of parton-medium inter- 
action and that pQCD based models are quite able to ac- 
count for the data. Given the lack of systematic studies 
with realistic medium geometries and the uncertainties 
still associated with a correct determination of parton 
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energy given a jet in a medium background, the observ- 
able has not been included into the analysis presented 
here. 

There is currently no reason to expect that for in- 
stance YaJEM-DE as constrained from Raa and Iaa 
would not be in agreement with the dijet asymmetry 
data. Qualitatively, the YaJEM-DE is very similar the 
scenario described in , i.e. medium-induced soft gluon 
radiation is depleted quickly by elastic processes whereas 
hard gluon radiation is more robust and survives as part 
of the jet. From the results of [ll| it is expected that 
medium-modified jets largely resemble jets in vacuum 
but shifted in energy. However, a (numerically rather 
involved) quantitatve comparison is still missing. 



E. Color (de-)coherence in the medium 

An essential ingredient to the detailed description of a 
developing parton shower in vacuum is color coherence, 
leading effectively to the phenomenon of angular order- 
ing, i.e. the angles between parent and daughter parton 
in subsequent gluon emissions decrease. For in-medium 
radiation on the other hand, i t ha s been found that anti- 
angular ordering is obtained jll6l | and that decoherence 
of the total rad iatio n spectrum in a sufficiently dense 
medium follows |ll7{ . 

Such a density-driven transition from color coherence 
to decoherence has so far not been included into any 
model of parton-medium interaction, thus all models cur- 
rently being compared to data miss a potentially impor- 
tant contribution to the dynamics. 

There are however two reasons to believe that this is 
not a huge problem in practice. First, if one estimates the 
formation time by the uncertainty relation as paramet- 
rically r ~ E/Q^ where E is the parton energy and 
the current virtuality scale, then the first hard branchings 
(which largely determine the structure of a jet in terms 
of e.g. subjet structure or transverse broadening) hap- 
pen at timescales much smaller than any medium forma- 
tion timescale as discussed in e.g. hydrodynamical initial 
states. This implies that whatever a dense medium will 
do to color coherence afterwards, the basic structure of 
a jet is determined by color coherence just as in vacuum 
since the medium cannot yet have formed. If a medium 
could be prepared beforehand and a jet embetted into 
this medium, then color decoherence would be a much 
more severe issue, but this is a different problem which 
does not correspond to the experimental situation. 

Second, while the transition to color decoherence is a 
difficult problem, the difference between full color coher- 
ence and full decoherence can be explored in in-medium 
MC shower models. This has been done e.g. in [t^ and 
was found to be a small effect as compared to other model 
uncertainties. Thus, while color decoherence is a very 
interesting theoretical problem which is not sufficiently 
addressed in current modelling and highly relevant for 
precision jet physics, there is currently no reason to as- 



sume that it would be a major effect in practice or that 
it would invalidate any of the results given above. 



F. Core-corona picture 

In order to explain the large 5*" ^ ob served at RHIC 
kinematics, core-corona models like jllSl ] have been pos- 
tulated in which the parton-medium interaction is essen- 
tially geometry driven: If, after a certain formation time 
Tforrm the parton is still in the interaction region, the 
parton is absorbed by the medium, otherwise the parton 
escapes without substantial interaction. 

In the context of angular dependent observablcs, such 
geometrical models work in general better than pQCD 
energy loss computations. However, they fail for other 
observables. In essence, a formation time in the energy 
loss picture can be written as a pathlength L dependence 

— with n a large number. Such a 

model would thus be disfavoured by the LHC data in the 
same way as the AdS dependence of energy loss. i.e. a 
purely geometric model could not generate sufficient Pt 
dependence to account for the data. 

Another problem is the act ual p arameter Tform needed 
to account for the data. In [ll8l | this is found as ~ 2.3 
fm, i.e. significantly longer than the formation times 
needed by typical hydrodynamical bulk models. While 
it is pointed out correctly that the hydrodynamical for- 
mation time Ti measures the time by which the medium is 
sufhciently equilibrated and pressure gradients act, and 
that this doesn't need to be equal to the time at which 
parton-medium interaction is strong, the difference un- 
fortunately acts in the wrong way: Parton-medium in- 
teractions probe the density of color charges, and color 
charges must be present in the system before it can equi- 
librate, thus Tform < Ti is the expected ordering. 

The observed form of /^^(^r) as discussed in llVDI like- 
wise argues against a geometrical picture in the sense that 
the enhancement seen at low zt is qualitatively and quan- 
titatively consistent with hadronizing medium-induced 
radiation. In a core-corona picture where a parton is ei- 
ther absorbed or escapes unmodified, no reason for such 
an enhancement would be present. 

Thus, while a geometrical interpretation is tempting 
(and indeed superior) given a subset of the high Pt data, 
it is insufficient to account even qualitatively for other 
parts of the data and has for this reason not been con- 
sidered in the main analysis section of this work. 



VI. CONCLUSIONS 

A systematic study of combinations of parton-medium 
interaction models and bulk medium fluid dynamical evo- 
lution models against a large body of available precision 
data on high Pt observables allows to establish the fol- 
lowing points: 
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• The basic strategy works, i.e. the multi-observable 
analysis has constraining power for both parton- 
medium interaction model and bulk medium evolu- 
tion model. The constraints posed by the data are 
significant, more than 95% of the scenarios tested 
showed a failure to describe the data for at least one 
observable. However, the strategy relies strongly 
also on the use of hydrodynamical modelling which 
in itself must be constrained against bulk observ- 
ables, it would fail if also parametrically plausible 
medium model choices without explicit constraints 
from the data (Bjorken cylinder, hard sphere over- 
lap...) would be acceptable. 

• The constraints on the parton-medium interaction 
side identify a number of features of the relevant 
physics process: The dominant energy transfer 
mechanism from the hard parton to the medium 
must be coherent. It leads to induced radiation 
which is observed in experiment and is qualitatively 
in agreement with medium induced gluon radiation 
in a pQCD scenario. A subdominant, elastic com- 
ponent is needed to deplete and decorrelate the soft 
gluon spectrum from the jet structure via energy 
and momentum transfer into the bulk medium, this 
component can be constrained from above and from 
below to about 10%. Presumably, there is an ex- 
plicit dependence on the initial hard parton energy 
in the MMFF, leading to P{AE,E) rather than 
P{AE) in the energy loss approximation. It follows 
from the above that modelling finite kinematics and 
full in-medium shower development is crucial to get 
agreement with all features of the data. 

• Constraints for the hydrodynamical side favour 
generally a late onset of energy loss and a large 
extension of the freeze-out hypersurface in space. 
Both a steeper initial density profile as found in 
CGC scenarios and viscous entropy generation are 
favoured by the data, but somewhat surprisingly 
are not dominant effects. In spite of commonly 
made assumptions, there is no reason to identify the 
equilibration time of the medium tq or the freeze- 
out isothermal surface characterized by Tp as the 
relevant scales for the parton-medium interaction. 

• Other high Pt observablcs so far do not substan- 
tially challenge this picture: Both heavy quark 
suppression and identified proton suppression are 
known to contain different and non-perturbative 
physics, as the hadronization cannot be argued to 



take place outside the medium, thus it is likely that 
the apparent puzzles associated with these observ- 
ablcs translate into constraints for this new ingredi- 
ent. Jet observables at the LHC are qualitatively in 
agreement with the picture outlined above, quan- 
titative comparisons are difficult and numerically 
challenging due to various jet-specific issues, but in 
many cases remain to be done. 
The weakest point of almost all models studied here 
remain observables which are differential in the reaction 
plane angle 0. It is in general non-trivial to reproduce 
the observed magnitude of S'™^ — even the calculations 
which are tuned to maximize this quantity barely reach 
above the measured value. Most probably this has to be 
be attributed to a problem on the hydrodynamical side 
rather than the parton-medium interaction side, since 
Iaa for instance also probes pathlength dependence of 
the medium effect, but shows no problem comparable to 
5™ J. This needs further study including initial state fluc- 
tuations of the bulk medium and different choices for the 
exit criterion of a parton from the medium. 

Since the strategy used in this work exploits qualita- 
tive differences to discriminate between various models, 
at this point no conclusion with regard to numerical val- 
ues of medium properties such as the transport coefficient 
q can be drawn. This is largely due to what has been 
summarized before as 'implementation details' and dis- 
cussed in depth in [t^: Internal model parameters such 
as cutoffs lead to a large uncertainty in the extraction 
of numerical values. It is unlikely that all this uncer- 
tainty can be overcome by a data driven approach, thus 
improved theoretical calculations based on models which 
contain the right qualitative dynamics is needed to turn 
hard probes into a quantitative measurement of medium 
properties. 
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